Thouless Energy and Correlations of QCD Dirac Eigenvalues 
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Eigenvalues and eigenfunctions of the QCD Dirac operator are studied for an instanton liquid 
partition function. We find that for energy differences SE below an energy scale E c , identified as 
the Thouless energy, the eigenvalue correlations are given by Random Matrix Theory. The value of 
E c shows a weak volume dependence for eigenvalues near zero and is consistent with a scaling of 
E c ~ l/L 2 in the bulk of the spectrum in agreement with estimates from chiral perturbation theory, 
that E c /A « F%L 2 jir (with average level spacing A). For SE > E c the number variance shows a 
linear dependence. For the wave functions we find a small nonzero multifractality index. 
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Random matrix theories have been applied to many 
aspects of mesoscopic systems (see for recent re- 

views). In particular, eigenvalue correlations have re- 
ceived a great deal of attention in this context. Two im- 
portant energy scales have been identified: i) the Thou- 
less energy, defined as the inverse diffusion time of an 
electron through the sample, i.e. 



HD 



(i) 



where D is the diffusion constant, and ii) the energy 
h/r e , where r e is the elastic collision time. Eigenvalue 
correlations on a scale SE can then be classified ac- 
cording to following three regimes: the ergodic regime 
for SE < E c , the diffusive or Altshuler-Shklovskii ||] 
regime for E c < SE < h/r e and the ballistic regime for 
SE > h/r e (see @j for recent work on this topic). 

The eigenvalue correlations can be measured conve- 
niently by the number variance, ^(n). This statistic is 
defined as the variance of the number of levels in an in- 
terval that contains n levels on average. In the ergodic 
regime, eigenvalue correlations are given by the invariant 
random matrix ensembles with ^(n) ~ (2/f3ir 2 ) log(n). 
In the diffusive regime the situation is more complicated. 
For a critical system, with a localization length that 
scales with the size of the sample, it is argued W that 
T,2(n) = x n (with x < I), whereas for weaker disorder 
the expectation is that £2(71) ~ n d / 2 . For a critical sys- 
tem, the slope of the number variance has been related 
to the multifractality index of the wavefunctions jq] . 

In this letter we wish to investigate to what extent 
such scenarios are realized in QCD. We will investigate 
eigenvalues of the Euclidean Dirac operator. Because 
of the J7a(1) symmetry they occur in pairs ±Afc or are 
zero. What is of main interest are the eigenvalues near 
zero which, for broken chiral symmetry, are spaced as 
A = tt/YiV (the space time volume is denoted by V). 
This is based on the Banks-Casher formula m according 
to which the order parameter of the chiral phase tran- 
sition, E, and the spectral density near zero are related 
by E = np(0)/V. It is therefore natural to define the 



microscopic limit, V — ► 00 with u = AVE kept fixed and 
the associated microscopic spectral density H 



Ps(u) = lim (p( — —)). 



(2) 



There is ample evidence from lattice QCD |9|JI0|] and in- 
stanton liquid simulations that ps(u) and other cor- 
relators on the scale of individual level spacings [|l2| are 
given by chiral Random Matrix Theory (chRMT), i.e. 
RMT's with the chiral symmetries of the QCD partition 
function. However, at scales beyond a few eigenvalue 
spacings in both instanton [jllj and lattice QCD [p||lC|] 
simulations the Dirac eigenvalues near zero show stronger 
fluctuations than in chRMT. This indicates the presence 
of an energy scale in QCD which may be identified as the 
Thouless energy. The interpretation of spontaneous chi- 
ral symmetry breaking as a delocalization transition was 
made earlier in [fl3f . By analogy with the Kubo formula, 
E plays the role of the conductivity |r| . 

Because of the chiral symmetry of the Dirac operator 
and its spontaneous breaking, the eigenvalue correlations 
near zero in the ergodic domain are given by the chiral 
ensembles [^,0- This is the domain where pion loops 
can be ignored. Its boundary is thus given by [ |l5| , |l6| a 
valence quark mass scale m c where the mass of the asso- 
ciated Goldstone boson, y/m c B (according to the PCAC 
relation B = Y./F 2 , where is the pion decay constant), 
is of the order of the inverse linear dimension of the box. 
This relation can be rewritten as flTnl 
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(3) 



It is therefore tempting to interpret 1/B as the diffu- 
sion constant. Because B is large on a hadronic scale 
(B w 1660 MeV), the diffusion constant is relatively 
small. With eigenvalue spacing A = 7r/EV, this con- 
dition can be rewritten in dimensionless form as 



mc = Fl L 2 

A 7T 



(4) 
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Here, g c plays the role of the dimensionless conductivity. 
Another discussion of F n in terms of the conductivity is 

for an Na 4 lattice, this 



given in [18 



relation reads g c — F 2 a 2 ^/N/n. On a 16 4 lattice with a 



In lattice QCD, 

w 

lattice spacing of 0.2 fm this results in a dimensionless 
Thouless energy of about one level spacing. The above 
discussion was for sea quark masses much less than the 
valence quark mass scale m c . Further subtleties arise in 
the quenched limit (see Q for more details). 

The chiral Random Matrix Theories for Nf massless 
quarks in the sector of topological charge v are defined 
by the partition function Sib! 



J N f M 



= I DWdet Nf 



iW 
iW^ 



-nf3TrV(W^W) 



(5) 



where W is a n x (n + v) matrix. The parameter 2n is 
identified as the dimensionless volume of space time. In 
this paper we only consider the chiral Gaussian Unitary 
Ensemble (chGUE) with complex matrix elements W ( 
(3 = 2) and a Gaussian potential V(x) = Y. 2 x. In this 
case ps(u) and all other correlation functions have been 
derived analytically (TJj2^,|l^] . We only note that in the 
bulk of the spectrum the correlations are given by the in- 
variant RMT's. It can be shown JUII that > for & = 2 > 
ps(u) and other correlators do not depend on the poten- 
tial V{x). The application of chRMT to QCD has been 
put on a firm foundation by these and other universality 
proofs p^-|2^]. Whether or not QCD is in this universal- 
ity class is a dynamical question that can only be proven 
by explicit numerical simulations. 

Our calculations will be performed using the instanton 
liquid model. In this model the gauge field configurations 
are given by a superposition of instantons. The Euclidean 
QCD partition function is then approximated by 



Dndet N f(D + m)e 



-5Vm 



(6) 



where the integral is over the collective coordinates of 
the instantons and the Dirac operator is denoted by D. 
For each instanton we have 12 collective coordinates (for 
three colors). The Yang-Mills action is denoted by Sym- 
The fermion determinant is evaluated in the space of the 
fcrmionic zero modes of the instantons. In our calcula- 
tions we use the standard instanton density N/V = 1 (in 
units of /m _1 ). For further discussion of this partition 
function, which obeys the flavor and chiral symmetries 
of the QCD partition function, we refer to pqj. 

The partition function (^) is evaluated by means of 
a Metropolis algorithm. We perform on the order of 
10,000 sweeps for each set of parameters. The eigenval- 
ues and eigenvectors of the Dirac operator are calculated 
by means of standard diagonalization procedures. In this 
letter we restrict ourselves to Nf = 0. This is a natural 
choice because our results are compared with ideas from 



disordered mesoscopic systems where no fermion deter- 
minant is present, and, moreover, it allows us to study 
much larger volumes. 

In order to separate the fluctuations of the eigenvalues 
from the average spectral density, the spectrum is un- 
folded such that the unfolded spectrum has unit average 
level spacing. All our spectral observables are calculated 
with the unfolded spectrum. 



■ N = 128 
• N = 512 
- chGUE 
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FIG. 1. The number variance £2(71) versus n approxima- 
tion for an interval starting at A = 0. The total number of 
instantons is denoted by N . 

In Figs. 1 and 2 we show the number variance, £2(72) 
versus n for various total numbers of instantons with in- 
stanton density N/V = 1. The chGUE result for S 2 (n) 
|l0f is depicted by the solid curve. In both figures, the 
upper figure (with only two volumes) is a blown-up ver- 
sion of the lower figure. In Fig. 1, T<2(n) is calculated for 
the interval starting at A = 0. Fig. 2 represents the num- 
ber variance in the bulk of the spectrum obtained from 
an interval that is symmetric about the average positive 
unfolded eigenvalue. In both cases we observe a clear 
transition point n c below which the number variance is 
given by RMT. In Fig. 1, the value of the crossover point, 
n c « 2, depends only weakly on the total number of in- 
stantons (or the volume). This is not in agreement with 
the theoretical expectation (Q) that n c sw F 2 W /it in 
four dimensions. For N/V = 1 fm~ A the value of n c for 
N instantons is given by n c w 0.07\/]V which is on the 
order of the results found in Fig. 1. However in the bulk 
of the spectrum (Fig. 2), the value of n c is consistent 
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FIG. 2. The number variance Y>2(n) versus n in the bulk 
of the spectrum. 



inverse participation ratio that is not too different from 
the random matrix result (full line). The eigenfunctions 
corresponding to small and large eigenvalues appear to 
be somewhat more localized. A more definitive result 
for the character of the wavefunctions follows from the 
scaling behavior of hW with the volume. A double log- 
arithmic plot of ATJ 2 (A) versus N is shown in Fig. 3 
(lower). Results are given both for the energy intervals 
[0.1,0.2] (open circles) and [0.7,0.8] (full circles). The 
first region corresponds to a part of the spectrum where 
the number variance shows a linear behavior for the vol- 
umes shown in Fig. 1 , and the second region corresponds 
to the bulk of the spectrum. The multifractality index r\ 
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with a y/V scaling but the numerical constant appears 
to be larger than the above estimate. This result is in 
agreement with the finding that correlations of lattice 
QCD eigenvalues are given by RMT up to distances of 
more than 100 spacings [|l2|j27| . 

Beyond the crossover point the number variance shows 
a linear behavior with a slope x ~ 0.08 for eigenvalues 
near zero and x ~ 0-04 in the bulk of the spectrum. The 
downward trend of the curves for larger values of n is a 
well understood finite size effect. This prevents us from 
saying more about the ballistic regime, an energy scale 
of roughly the inverse distance between instantons. 

In the ergodic regime we expect that eigenvalue corre- 
lations are given by the chiral random matrix ensembles. 
Indeed, both the microscopic spectral density up to two 
level spacings and the nearest neighbor spacing distribu- 
tion are in perfect agreement with the chGUE. 

The number of significant components of the wave 
function is measured by the participation ratio. Its in- 
verse is defined as 



(7) 



The tpkW are the normalized A-component eigenfunc- 
tions (corresponding to eigenvalue A) of the Dirac op- 
erator in the space of the fermionic zero modes of the 
individual instantons. Results for A/2 (A) as a function 
of A are depicted in Fig. 3 (upper). The general im- 
pression is that the wave- functions are extended with an 
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FIG. 3. The inverse participation ratio times N versus the 
corresponding Dirac eigenvalues (upper), and the scaling be- 
havior versus N (lower) for two energy intervals. 



is defined by 



(8) 



According to H the value of rj = 2x,d (where x is the slope 
of the linear piece in the number variance) in the criti- 
cal domain. From the volume dependence of I2 shown in 
Fig. 3 (lower) we find values for r\jd of about 0.02 and 
0.04 for the intervals [0.1,0.2] and [0.7,0.8], respectively. 
These values are well below the theoretical result of 2%. 
Apparently, our ensemble of instantons is not in the criti- 
cal region. Our results should be contrasted with Wilson 
lattice QCD Dirac eigenfunctions which were found to be 
localized p8| . We have no explanation for this discrep- 
ancy. 
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In conclusion, we have identified an energy scale be- 
low which the eigenvalue correlations of the QCD Dirac 
operator are given by chRMT. In analogy with the the- 
ory of mesoscopic systems, this scale will be called the 
Thouless energy. For eigenvalues near zero we find a 
Thouless energy that only shows a weak volume depen- 
dence, whereas for eigenvalues in the bulk the Thouless 
energy scales roughly with the square root of the volume 
in agreement with theoretical prediction. 

For energy scales beyond the Thouless energy a lin- 
ear n-dependence of the number variance is found. Our 
wave functions show a small nonzero multifractality in- 
dex which does not obey the relation derived for crit- 
ical mesoscopic systems. Interesting connections with 
the scalar susceptibility and quenched chiral perturba- 
tion theory will be discussed elsewhere [^9| . 

A. Smilga and Y. Fyodorov are thanked for stimulat- 
ing discussions. After the completion of this work we 
received a preprint (R. Janik et al., bep-ph/9803289 ) in 
which similar ideas were discussed. In that work the 
diffusion constant is related to diffusion in a 4+1 dimen- 
sional space time. 
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